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Abstract 

This paper deals with the effect of shaft frequency on the performance of 
a submerged journal undergoing noncentered circular whirl. The main emphasis 
of the paper Is on the behavior of the vapor cavitation bubble and Its effect 
on the bearing performance as a function of frequency. A cavitation algorithm 
due to Elrod was Implemented In a computer program which solves a 
time-dependent Reynolds equation. This algorithm automatically handles the 
boundary conditions by using a switch function and a control volume approach 
which conserves mass throughout the entire flow. The shaft frequencies In 
this Investigation ranged from 0 rad/s (squeeze-film damper) to -104 rad/s (a 
case In which oil-whip condition was produced momentarily). For the 
particular vibration amplitude chosen In this Investigation It was observed 
that vapor cavitation had an effect on the load components for the full range 
of shaft frequencies Investigated. 

Nomenclature 

D shaft diameter , m 

e eccentricity, m 

e^ dynamic eccentricity, m 

e $ static eccentricity, m 


f R * racial load component, N 
F tangential load component, N 

<p 

g switch function (cavitation Index) 

h film thickness, m 

L/D length to diameter ratio 

M bearing torque, N-m 

m lineal mass flux, kg/m-s 

0^ bearing center 

journal center 

2 

p fluid pressure, N/m 

2 

p a ambient pressure, N/m 

a 

2 

P c cavitation pressure, N/m 

R radius of shaft, m 

t time, s 

U sum of the surface velocities In x-dlrectlon, m/s 

Y sum of the surface velocity vectors, m/s 

W load capacity, N 

w squeeze-velocity, m/s 

2 

(1 liquid bulk modulus, N/m 

Y angular position of minimum-film thickness, rad 

Y angular velocity of the minimum-film position, rad/s 

AR radial clearance, m 

c total eccentricity ratio, e/AR 

e d dynamic eccentricity ratio 

c $ static eccentricity ratio 

p fluid density In the fluid film, kg/m 3 

p c fluid density within cavltated zone, kg/m 
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e fractional film content in cavitation zone, density ratio ( P / Pc ) in 

full-film zone 

p dynamic viscosity, N-s/m 

period for 1 revolution of the journal about its own center, 2ir/u s ,s 

Td period for 1 revolution of the journal center about a fixed point 
relative to housing center, 2ir/w ( j, s 

(p angular coordinate relative to minimum-film line, rad 

«1» attitude angle, tan -1 (F /F D ), rad 

cp K 

wj orbital angular velocity of journal center about a fixed point relative 
to housing center, rad/s 

cj $ angular velocity of journal about its own center, rad/s 
INTRODUCTION 

This paper deals with the effect of shaft frequency on the performance of 
a dynamically loaded journal bearing with vapor cavitation. Despite the 
Importance of cavitation in dynamically loaded bearings, very little work has 
been reported which fully accounts for the formation of vapor cavitation. 
Because of the complexity and the computational requirements of the process, 
most researchers are compelled to assume that in a dynamically loaded bearing 
the boundary conditions remain similar to a steadily loaded journal bearing 
(1). An earlier paper by Brewe (2) Illustrated some of the consequences of 
not taking film reformation and the moving boundary into account. 

Cavitation can manifest Itself in two distinct forms: gaseous cavitation 

and vaporous cavitation. Gaseous cavitation is generally seen in a steadily 
loaded journal bearing. It occurs in low pressure regions of the flow where 
the saturation pressure of the fluid has been reached. This form of 
cavitation is benign to bearing surface damage (3). In contrast, vaporous 
cavitation normally occurs in dynamically loaded bearings and is responsible 
for physical damage to the bearing surfaces as a result of the collapse of the 
vaporous cavities on the surfaces. Jacobson and Hamrock (4,5) recorded the 
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appearance of vapor cavitation In a dynamically loaded journal bearing for the 
first time by using high-speed photography. More recently, Walton et al. have 
used high-speed photography to observe cavitation In squeeze-film dampers. 
("Experimental Observation of Cavltatlng Squeeze Film Dampers" presented at 
the 1986 ASME/ASLE Tribology Conference, Oct. 20-22, 1986, Pittsburgh, PA. To 
be published In Journal of Tribology, ASME.) 

Further, cavitation plays an Important role In the stability of rotating 
machinery. A previous investigation by Brewe (2) verified that cavitation can 
generate some load capacity and should not be taken lightly in theoretical 
studies. Therefore, It Is Important to further study the effect of cavitation 
on the performance of the bearing. 

This study centers on the theoretical investigation of the performance of 
the bearing with particular attention to cavitation effects. The cavitation 
algorithm used Is based on the mass conservation throughout the entire flow 
region. 

THEORY 

Cavitation Algorithm 

A computer program was developed based on Elrod's algorithm (6) that 
combines a control volume approach to deriving the finite-difference equation 
with a switching function that automatically handles the moving boundary 
condition at the liquid-vapor Interface. The complication In satisfying the 
boundary condition Is avoided by defining a single Independent parameter e 
which links the full-film region and the cavltated region. This parameter is 
defined by the ratio of the local density to that of the cavitation density, 

1 .e. , 

e = *- 
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It follows that In the cavitation region e < 1.0 and in the full-film 
region e > 1.0. A switch function g was defined as "on" (g = 1) and 
"off" (g = 0) In the full-film and cavitation regions respectively. By using 
this property, a universal Reynolds equation was defined as Indicated: 

^t^ + * * m = 0 


where m 


represents the lineal mass flux and Is given by 


-v 

m 


ehv 

2 



(ig(e)ve 


[ 2 ] 


In the cavitation region, the second term vanishes since g = 0. This 
equation was Integrated In time by using an Euler's method to advance in 
time. At every time step the difference equation was solved by an alternating 
direction Implicit (ADI) scheme. The switch function was updated after each 
half-time step. Occasionally, during the collapse of the vapor bubble, liquid 
was piling up at the boundary. This was due to an Inherent time lag of the 
switch function with the calculation. Consequently, the switch function was 
Immediately updated to accommodate the movement of the boundary, and the 
calculation was reiterated. 


The results of the program developed were authenticated by comparison to 
an experiment of Ref. (4). When compared to the pseudo-Gumbel boundary 
conditions (2), the results predicted a significant difference in the computed 
load capacity (up to 20 percent). Furthermore, It was found that the size, 
life, and location of the cavitation bubble were considerably affected by 
imposing the pseudo-Gumbel boundary conditions. 
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Fluid Velocity Under Dynamic Loading 


Consider the bearing geometry shown In Fig. 1. As shown, the minimum-film 
thickness rotates clockwise around the bearing at a whirl frequency of y. 
Jacobson and Flamrock (5) show that from this diagram y can be evaluated 

c. sln(w.t) 

f 3 1 

+ c d cos(w d t) 1 

so that 



C{j [c s COS(w d t) + c d ] co d 
Y = ~2 2 

c + g. + 2 c e . COS( W . t ) 
s d s d d 

The total eccentricity e defined as the Instantaneous distance between 

0. and 0, Is: 
b j 


e = 


e d s1n < w d t} 


sin y 

The film thickness Is, therefore, defined as follows: 

h = AR( 1 - c cos tp) 


[5] 

[ 6 ] 


The sliding velocity of the journal Is 


U = (u $ - 2 y) R 


[7] 


Test Bearing 

The bearing simulated In this study Is based on an experimental study of 
a journal bearing submerged In Dexron lubricant (4). The journal rotated 
about Its own fixed center while the bearing vibrated In a noncentered 
circular motion keeping the axes of vibration parallel to that of rotation. 

Note that there are two time scales In this problem: one, t , Is associated 
with the spin of the journal about Its own center; the other, x d , Is associated 
with the translational motion of the journal center about a fixed point. The 
ratio of these time scales T ( j/ T s is, In some sense, a measure of wedge to 
squeeze action taking place. This range varied from 0.00 to 1.13 In this 
study. The data used In the simulation are summarized In Table 1. 
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RESULTS AND DISCUSSION 

Effect of Shalt Frequency on Cavitation 

Table 2 lists the various shaft frequencies run In this analysis and 
Indicates the effect they had on cavitation. The times listed in the table 
are relative to the start of the orbit. The period for the noncentered shaft 
whirl was 67.78 ms for each case. Cavitation persisted for a longer period of 
time in the case of the squeeze-film damper (i.e., 48.91 ms) than in the other 
cases listed. The bubble life decreased with increasing shaft frequency. The 
cavitation desinence (bubble collapse) was weakly affected by shaft frequency. 
However, as the shaft frequency was Increased the cavitation inception was 
delayed which gave rise to a shorter bubble life. 

Orbital Pressure and Cavitation History 

Figures 2 and 3 Illustrate the responding pressure variation to the 
motion of the shaft undergoing noncentered circular whirl for two contrasting 
situations. The case of a squeeze film damper (i.e. = 0) in Fig. 2 Is 

contrasted with the case in which the oil-whip condition Is momentarily 
observed (i.e. o> s = -104.29 rad/s) in Fig. 3. Both cases have an orbital 
frequency of -92.7 rad/s. In each figure, part (a) illustrates the starting 
position (i.e. e = 0.1). Parts (a) to (f) depict several key events in the 
clockwise circular motion of the journal center within the bearing. 

Figure 2(a) represents the position of the journal within the housing and 
the associated pressure distribution at the instant that the eccentricity is 
at a minimum. The dark oblong region indicates the shape and position of a 
vapor bubble downstream of the minimum-film line. The position of the bubble 
within the small clearance space of the journal bearing is also indicated. At 
this particular instant the tendency of the pressure flow along the side of the 
bubble to cause the bubble to collapse becomes greater since it is proportional 

3 

to h . The journal center moves in a clockwise orbit from Figs. 2(a) to (b) 
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(nearly one-eighth of a revolution). While the eccentricity Is Increasing, 
the minimum film Is decreasing which gives rise to the squeeze action In the 
vicinity of the minimum-film line. A slight rise In the pressure field can be 
observed there. Note In Fig. 2(b) that the vapor bubble has collapsed farther 
downstream from Its position In Fig. 2(a). It Is located closer to the 
position of maximum-film thickness. 

Cavitation Is absent from Fig. 2(b) to (c). In Fig. 2(c), a new 
cavitation bubble appears, and the pressure generation in the region of the 
minimum-film line becomes significantly large because of the squeeze and wedge 
effects taking place there. Figure 2(d) represents the configuration of the 
bearing when the maximum pressure occurs. Note that this happens 5.3 ms 
before the maximum eccentricity is reached because of the motion of the shaft 
[Fig. 4(a)]. As the journal approaches the end of the squeeze cycle and begins 
to separate from the shaft, the sliding velocity approaches a minimum value and 
the squeeze velocity passes through zero. The effect of this on the pressures 
more than offsets any increase In pressure that would have been realized 
because of the smaller minimum-film thickness. 

Figure 2(e) represents the end of the squeeze cycle and the beginning of 
the separation cycle. Now the minimum film begins to increase which creates a 
suction effect, an accompanying growth in the bubble, and a dissipation of the 
pressure hump noted In Fig. 2(f). It is interesting to note that the 
eccentricity is the same (0.78) for Figs. 2(d) and (f). The drastic difference 
in the pressure plot is a result of the squeeze velocity having opposite signs 
in the two situations. The minimum-film thickness continues to increase as it 
moves in a clockwise direction. The representation shown in Fig. 2(a) is 
repeated as are the representations for each succeeding orbit. 

A very significant difference between Figs. 2(e) and 3(d) occurs when the 
eccentricity is a maximum (i.e. e = 0.80). In the case of the squeeze-film 
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damper (i.e. o> s = 0), Fig. 2(e) Illustrates a sizeable pressure buildup and 
a well developed cavitation bubble. However, for the case in which 

= -104.29 rad/s. Fig. 3(d) shows no pressure buildup or cavitation. This 
condition arises because the net flow produced by the squeezing and sliding 
motion Is zero. The whirl frequency y at this Instant Is equal to one-half 
the shaft frequency and Is commonly referred to as the "half-frequency whirl" 
condition. 

Load as a Function of Shaft Frequency 

Figure 5(a) to (d) Illustrates how the load- carrying capacity varies as a 
function of the eccentricity ratio for various shaft frequencies . A case 

representing a squeeze-film damper (u^ = 0) is shown In Fig. 5(a). The load 
variation through one orbit of circular whirl resembles a hysteresis effect 
which is due to the occurrence of cavitation In the cycle (7). This effect is 
prevalent for all of the shaft frequencies studied In this Investigation. 

Each portion of parts (a) and (d) of Fig. 5 is labeled for discussion 
purposes; the reader should simultaneously refer to the corresponding pressure 
distributions presented In Figs. 2 and 3, respectively. 

For the particular cases studied, it can be said that increasing the 
shaft frequency had little effect on the onset of deslnation (bubble collapse) 
and the location of the peak load (viz. Fig. 5). However, the magnitude of 
the peak load decreased with the increase In the shaft frequency. 

The variation of the shape of the W-e curve with frequency is 
Interesting and well worth discussing. The portion of the cycle beginning 
with the squeeze motion through the peak pressure [Fig. 5(a) to (d)] is 
essentially invariant with the shaft frequency. But around the vicinity of 
c = 0.8 where the separation occurs, the effect of shaft frequency is more 
apparent. Here the W-e develops an increasingly steep gradient with 
increasing shaft frequency. Finally for the case = -104.29 rad/s 
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[Fig. 5(d)] the load capacity vanishes completely when the whirl frequency (y) 
Is equal to one-half the shaft frequency. This frequency can be predicted by 
examining equation [7]. The squeeze and sliding velocity for this case Is 
plotted In Fig. 4(c) . 

It may also be noted that on the return cycle (Fig. 5 from point f to a) 
there Is no appreciable load variation where = 0. As the frequency 
Increases, the variation In load with eccentricity In that portion Increases. 
This Is a noticeable variation (roughly 200 N) at the shaft frequency of 
-104.29 rad/s. After passing through the load reduction, the sliding and 
squeeze term In the Reynolds equation recovers to contribute to the load 
capacity. Another contributing factor to the load capacity in the separation 
cycle Is the cavitation. This can occur when the bubble moves downstream from 
the minimum-film line so that It positions Itself diametrically opposite to 
the directional load line described by the attitude angle. 

Radial and Tangential Loads 

Figure 6 illustrates how the shaft frequency can affect the dynamic 
loading (radial and tangential) as a function of eccentricity through one 
whirl orbit. Both the radial and tangential loads demonstrate a "hysteretic- 
llke" response with eccentricity. This effect is a result of cavitation 
occurring through some portion of the cycle. [See Refs. (2) and (7).] 
Operation at the different shaft frequency conditions has a very pronounced 
effect on the tangential load but only a slight effect on the radial load. As 
the shaft frequency is Increased, the tangential load plots become less 
hysteretic (See Fig. 6.) in nature. This is mainly because the onset of 
cavitation has been delayed at the higher shaft frequencies. Further, at the 
end of the squeeze cycle (i.e. e = 0.8 in this study) the tangential load 
decreases as the shaft frequency is increased. This is because at this 
particular stage in the motion the conditions for half-frequency whirl are 
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being approached for increasing shaft frequency. When the conditions for 
half- frequency whirl are met [i.e. in Fig. 6(c)], the tangential load and 
radial load are identically zero. 

Attitude Angle 

The effect of the shaft frequency on radial and tangential load is 
summarized in the attitude angle. The variation of the attitude angle with 
eccentricity through one circular orbit is shown in Fig. 7. These plots 
correspond to the case of a squeeze-film damper (« s = 0) and the case of 
synchronous shaft and vibration frequency (I.e. respectively. 

In the damper case [Fig. 7(a)], the asymmetry in the plot is again a result of 
the cavitation present throughout a portion of the whirl cycle. In Fig. 7(b) 
the delay of the onset of cavitation diminishes the asymmetry of the plot. 
Also, the magnitude of the attitude angle remains small (i.e. q/ < 30°) 
throughout most of the cycle (I.e. 0.15 < e < 0.79). At the end of the 
squeeze cycle the sign is changed, and it very quickly becomes small again for 
the remainder of the whirl cycle. In this instance the effect of the higher 
shaft frequency has been to suppress the magnitude of the attitude angle 
through a large portion of the cycle. 

Bearing Torque 

The variation of the bearing torque with eccentricity is shown in Fig. 8 
for three different cases. The bearing torque is a measure of the force 
required to overcome the friction force in the fluid film [See ref. (2)]. For 
the squeeze-film damper [Fig. 8(a)], the bearing torque is generally higher 
throughout the whirl orbit than in any of the cases that included shaft spin. 
As the shaft frequency Increased, the resultant effective sliding velocity 
decreased (Fig. 4) and, consequently, the amount of shear in the fluid. Note 
that the bearing torque becomes zero for the oil-whip condition [Fig. 8(c); 

C = 0.8]. 
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CONCLUSIONS 


A study of the effect of shaft frequency on bearing and damper 
performance under dynamic conditions Is presented. The motion considered In 
this study was that of noncentered circular whirl. The study revealed the 
following: 

1 . Increasing the shaft frequency caused a delay in the onset of 
cavitation. Since the collapse of the bubble occurs at nearly the same time 
in each case, the life of the bubble is decreased with increasing shaft 
frequency. 

2. During separation (i.e., 3h/3t > 0) Increasing shaft frequency 
results In an appreciable load variation with decreasing eccentricity. This 
is attributable to the size and location of the cavitation bubble and to the 
recovery of the wedge and squeeze effects. 

3. Operation at the different shaft frequencies has a very pronounced 
effect on tangential load, but only a slight effect on radial load. 

4. At the oil-whip condition it was verified that the bearing lost all 
load-carrying capacity, and bearing torque was zero. 
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TABLE 1. - OPERATING CONDITIONS 


AR, m 5.0 • 10 -4 

R, m 0.0425 

L/D 1.0 

E 0.1 to 0.8 

<■>(!• rad/s -92.7 

<o s > rad/s 0 to -104.29 

6 , N / m 2 1.72 - 1 0 9 

v, N-s/m 2 0.066 

P a . N/m 2 1.0133 • 105 

P c . N/m 2 0.0 


TABLE 2. - EFFECT OF SHAFT FREQUENCY ON CAVITATION 


Shaft 

frequency, 

rad/s 

Cavitation 

inception 

Cavitation 

desinence 

Bubble 

life, 

ms 

Time, 

ms 

c 

Time, 

ms 

c 

0.00 

24.300 

0.724 

73.210 

0.222 

48.91 

- 15.45 

26.205 

.751 

73.410 

.228 

47.21 

- 30.90 

28.305 

.774 

73.610 

.234 

45.31 

- 46.35 

30.305 

.789 

73.810 

.241 

43.51 

- 61 .80 

32.005 

.797 

74.010 

.248 

42.01 

- 77.25 

33.505 

.800 

74.110 

.250 

40.61 

- 92.70 

34.505 

.800 

74.210 

.254 

39.71 

- 104.29 

34.805 

.800 

73.910 

.243 

39.11 
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F I Lfi FILM FILM 


(a) BEGINNING OF SQUEEZE CYCLE. (b) BUBBLE COLLAPSE, (t = 6.2 MS; (C) PEAK PRESSURE, (t = 26.31 MS; 

(t - 0 ms; £ = 0.1.) £ = 0.24.) £ = 0.75.) 



FILM FILM FILM 

(f) BUBBLE EXPANSION, (t = 41.41 ms; (e) CAVITATION INCEPTION, (t = 34.81 MS; (d) BEGIN SEPARATION OIL-WHIP CONDITION. 

£ = 0.75.) £ = 0.799.) (t = 33.91 MS; £ = 0.80.) 

FIG. 3. - PRESSURE DISTRIBUTION AND BEARING CONFIGURATION FOR FULL PERIOD OF SHAFT WHIRL. (U g = -104.29 rad/s; FIGURES (a) TO (f) ARE 
VIEWED CLOCKWISE IN TIME.) 
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16. Abstract 

This paper deals with the effect of shaft frequency on the performance of a sub- 
merged journal undergoing noncentered circular whirl. The main emphasis of the 
paper Is on the behavior of the vapor cavitation bubble and Its effect on the 
bearing performance as a function of frequency. A cavitation algorithm due to 
Elrod was Implemented In a computer program which solves a time-dependent 
Reynolds equation. This algorithm automatically handles the boundary conditions 
by using a switch function and a control volume approach which conserves mass 
throughout the entire flow. The shaft frequencies In this Investigation ranged 
from 0 rad/s (squeeze-film damper) to -104 rad/s (a case In which oil-whip condi- 
tion was produced momentarily). For the particular vibration amplitude chosen In 
this Investigation It was observed that vapor cavitation had an effect on the 
load components for the full range of shaft frequencies Investigated. 
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